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Abstract. We consider the quantum Calogero model, which describes A'^ non- 
distinguishable quantum particles on the real line confined by a harmonic os- 
cillator potential and interacting via two-body interactions proportional to the 
inverse square of the inter-particlc distance. We elaborate a novel solution al- 
gorithm which allows us to obtain fully explicit formulas for its eigenfunctions, 
for arbitrary coupling parameter and particle number. We also show that our 
method applies, with minor changes, to all Calogero models associated with 
classical root systems. 



1. Introduction 

In this paper we elaborate a novel solution method for the iV-body Calogero 
model defined by the Hamiltonian 

N 

(1) ^ = E (-^'. + + 2A(A - 1) E 7^^^' 

where A > is the coupling parameter, Xj e R the particle coordinates, dx '■= 
d/dxj, and N = 1,2,3,... the particle number (we set the harmonic oscillator 
frequency w > to 1 without loss of generality: this parameter can easily be 
introduced by scaling xj \/uJXj, H uiH, etc.). As is well-known [Cal71,Sut72], 
this model has exact eigenfunctions of the form 

(2) Vn = V-oPn, 

where 

N 

(3) V'o(x)-ne~^'^'n(^^-^^)^ 

3=1 j<k 

is the groundstate eigenfunction and -Pn(x) are polynomials which are symmetric, 
i.e., invariant under permutations of the particle coordinates. These polynomials 
are labeled by A^-tuples n — {ni, . . . , njv) of non-negative integers, nj S Ng. Due 
to the permutation symmetry these labels can be restricted to partitions, i.e., 

n-i > n2 > . ■ ■ > TIN > 0, 

but we will not always make this restriction. The corresponding exact eigenvalues 
are given by the following remarkably simple formulas: 

(4) = 2(^1+712 + ••• + njv) + £;o, E^ = N{l + \{N -I)). 

We refer to the Pn as reduced polynomial eigenfunctions of the Calogero model, 
and our aim is to derive explicit formulas for them. These polynomials are a 
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natural many-variable generalization of the Hermite polynomials to which they 
reduce in the special case = 1 [BF97]. Previous results on these functions 
[Cal69,Sut72,Per71,Gam75,BHV92,Kak96,DLM04] will be discussed in more detail 
below. We also mention that there has been considerable interest in many-variable 
generalizations of classical orthogonal polynomials in the mathematics literature; 
see e.g. [DX01,Mac95] and references therein. 

Calogero found in his seminal paper [Cal71] the exact eigenvalues of a closely 
related model which differs from the one above only in its center of mass mo- 
tion (for the convenience of the reader we discuss the precise relation of these two 
models in Appendix . The eigenvalues of the Hamiltonian in |^ was given by 
Sutherland [Sut72], who also presented an algorithm for constructing the reduced 
eigenfunctions P^. This algorithm starts with the ansatz in which converts 
the eigenvalue problem for H into a problem of diagonalizing a certain triangular 
matrix. Thus, the eigenvalues of H can be read off from the diagonal of this matrix, 
and the eigenfunctions are determined by certain recursion relations which truncate 
after a finite number of steps; see Section 13 for details. To our knowledge, these 
recursion relations have not been solved by a closed formula. 

In Sutherland's paper [Sut72] the emphasis was on a translation invariant iV- 
body model with a 1/ sin^-interaction which is also exactly solvable and the solu- 
tion algorithm was elaborated in detail only for this so-called Sutherland model. 
In [LanOl] one of us presented an alternative algorithm to solve the Sutherland 
model which, different from Sutherland's, also can be generalized to the elliptic 
case; see [Lan05] and references therein. In the present paper we extend this solu- 
tion algorithm to the Calogero model. We stress that the Calogero model is more 
complicated than the Sutherland model due to the presence of the harmonic oscil- 
lator potential, and this leads to various interesting and novel features. It is also 
interesting to note that, in our approach, the factorization of the eigenfunctions in 
(01 is a consequence, rather than an essential ingredient, of the method. Moreover, 
rather than constructing the eigenfunctions as linear combinations of the free boson 
eigenstates as Sutherland, we obtain a set of somewhat more complicated functions 
which lead to simpler recursion relations which we solve explicitly. This gives our 
main result: an explicit formula for the reduced polynomial eigenfunctions of the 
Calogero model. 

We now briefly describe this result. For each fixed x e and e > 0, let Cj 
denote the following set of nested circles in the complex plane: 

(5) Cj : yj = ( max (Ixfcl) + eJ)e'^^ -tt < (p^- < tt, j = 1, . . . , N. 

l<k<J\ 

Using these curves as integration paths, define for each n G Nq the functions 




(6) /n(x) :^ M ( ^ I 



which are symmetric polynomials independent of e > 0; see Section |31 Our main 
result is a fully explicit formula for the functions Pn as linear superpositions of these 
functions /„• We use the natural basis elements e Nq defined by {ej)k ■= Sjk 
and write 5n(ni) := 5n,m for the Kronecker delta. 
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Theorem 1.1. For n e 'Z^ let 

m 

with the functions fm defined in and the coefficients 

oo ^ oc 

a„(m) = <5„(m) + E il^i E ' ' ' E E 

^y-j S = l jl<kl js<fcs l/l,...,!^s=0 

r=l 

where we use the shorthand notations 

(8) 5j7c(i^;m)= 2A(A - l)i^(l - 5j7c) - + l)5i/o^ifc, 

= nij + \{N + 

and 

(9) = {1 - ^)e, + {1 + iy)ek. 

Then Pn *s « reduced polynomial eigenfunction of the Calogero model corresponding 
to the eigenvalue in 

(The proof will be given in Sections 13. 11 - 13.41 ) 

It is important to note that the sums in Q only contain a finite number non-zero 
terms. It is also remarkable that Theorem ll.ll is non-trivial already for the simplest 
case TV = 1, as discussed in Section lXBl In Section 3.6 we use this result to construct 
somewhat more complicated basis functions than the /n, leading to another explicit 
formula for the reduced polynomial eigenfunctions; see Theorem 13. II 

Observe that for > 1, this result gives too many eigenfunctions: they are 
in Theorem 11.11 labeled by elements in Z^, but it is known that a complete set 
of eigenfunctions can be parameterized by partitions alone. Using the symbolic 
programming language MATHEMATICA we have checked for N — 2 that the 
are non-zero eigenfunctions also for non-partitions n, and we conjecture this to 
be true for all N. This over-completeness of our solution poses some interesting 
questions discussed in Remark 15. II 

Similarly as in Sutherland's algorithm [Sut72], we obtain the coefficients an(m) 
by diagonalizing a certain triangular matrix, and they are therefore non-zero only 
for m ^ n in some partial ordering ^; see Sectional However, this partial ordering 
is different from Sutherland's, and the matrix we get is simpler, which is why we 
can find its explicit eigenvectors. 

As found by Olshanetsky and Perelomov [OP77], the model discussed so-far can 
be naturally associated with the root system Ajv_i, and there are exactly solvable 
variants of the Calogero model related to all other root systems; see [OP83] for 
a comprehensive review. In particular, the Calogero models associated with the 
remaining classical root systems [OP83] can all be brought to the form of the _Bjv 
Hamiltonian 




(10) hb^^^i -e + ^; + "^"2 + 4A(A - 1) 
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with two coupling parameters /i, A > 0. To demonstrate the generahty of our 
sohition method we show that the construction of eigenfunctions with minor changes 
goes through also in this case, and we thereby obtain explicit formulas for a many- 
variable generalization of the Laguerre polynomials. This adds support to our hope 
that the method can be used so solve any Calogero-Sutherland type model. We 
should mention that the Bm Calogero model also can be solved using Sutherland's 
method, of course. 

As mentioned in the first paragraph, various other explicit results for the reduced 
polynomial eigenfunctions of the Calogero model exist in the literature. Calogero 
obtained such results for the cases TV = 2,3 [Cal69]. By exploiting an underlying 
group structure of the Hamiltonian, Perelemov [Per71] for iV = 4 and Gambardella 
for iV = 5 [Gam75] obtained the eigenfunctions in terms of "raising" operators 
acting on the groundstate. More recently these operator solutions were generalized 
to all N [BHV92, Kak96]. We also mention that Desrosiers et.al. obtained explicit 
results for the eigenfunctions of a supersymmetric generalization of the Calogero 
model using a determinantal construction [DLM04] . Our results seem different and 
complementary to these. 

The plan of the rest of this paper is as follows. In Section [21 we fix our notation 
and shortly review Sutherland's solution of the Calogero model [Sut72] and a simple 
variant thereof which, as we argue, is somewhat more natural. In Section |31 we 
present our solution of the A^-i Calogero model and thereby prove Theorem ll.il 
We also comment on the one-particle case, and we sketch a variant of our solution 
method which provides another explicit formula for the eigenfunctions. Our solution 
of the Sjv Calogero model is presented in Section^ We end with a few concluding 
remarks in Section |S1 Some technical details are deferred to two appendices. 



In this section we fix our notation and, to put our work into context, briefly 
review Sutherland's solution of the Calogero model [Sut72]. We will actually discuss 
a somewhat simpler variant of this solution method, as explained below. 

In the discussion below we make use of some notational conventions from the 
theory of partitions which we now recall; see e.g. [Mac95]. For partitions n = 
(rti,rt2, . . . ,n]\[), the non-zero are called the parts of n, and we use the short 
hand notation 



We also introduce a partial ordering of partitions: for two partitions m, n we write 



We will furthermore write m < n if m < n and |m| ^ |n|. 

The starting point of Sutherland's algorithm is the observation that the func- 
tion -00 in Q is the groundstate of the Hamiltonian in a fact which can 
be proved by a straightforward computation; see Remark 13. II As previously men- 
tioned, another key insight is that any eigenfunction of the Hamiltonian H in ^ 
can be factorized into a symmetric polynomial and the groundstate. This implies 



2. Sutherland's solution algorithm 



|n| := ni + 712 + ■ • • + ?^A^• 
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that such a symmetric polynomial is an eigenfunction of the differential operator 

N ^ 

(11) H := ^o'HiPo -Eo = Y. [-^X + 2^^^-.) " 2A5] ^^(9,^ - 9,J. 

3 = 1 j<k ^'^ 

The idea is now to construct these polynomials as linear combinations of the mono- 
mials 



(12) M„= ^P(i)---^ 



P{N) 



where n is a partition of length N and Sn the permutation group of N elements. 
We note in passing that the standard normalization of these monomials is different 
in that the sum in H12I) is restricted to the distinct permutations of the parts rij 
(see e.g. [Mac95]), but for our purposes the normalization where one sums over all 
permutations is more convenient. To proceed we use the fact that 

{-dl + 2x9^)0;" = 2na;" - n{n - l)a:"-2 

as well as the identity 



n—m—L 

— {d^ - dy){x''y"-' + 2/"a;™) = (n - m) ^ ^n^i-kyrn-i+k 

y fe=i 

-m{x"-'^y"'-^ + y"-ia;™-i), 



valid for all x,y € M. and to, n G No such that n > m. A proof of this identity can 
be found in Appendix B. It follows that 



(13) 



JV 

HMn = 2\n\Mn - ^nj(rij - l)M„-2e, 

— ^ (2 - (52iy,nj-nfc) ((1 — (5jy,o)n-j — n-fc) A^n~(i/+l)ej+(iy-l)efc 

where ['T-/2J = n/2 or (n — l)/2 for even or odd integers n, respectively. Hence, the 
action of H on the monomials Mn has triangular structure in the following sense: 

HMn = 2|n|M„ + bnmMm 

for certain coefficients &nm which can be determined from (|13|l . This suggests that 
H has eigenfunctions of the form 

(14) F„ = M„ + ^ MnmAfm 

m<n 

with corresponding eigenvalues — 2|n|. Indeed, inserting this result into the 
Schrodinger equation HP^ — EnPn and using the fact that the monomials Afm are 
linearly independent we obtain the following system of equations: 

(15) {En - Ejn)Unm = ^nm + ^ Wnk&km, HI < n. 

m<k:<n 
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It is important to note that |m| < |n| for all m < n. Moreover, for each partition 
n, there exists only a finite number of partitions m < n. Thus, (|15|l gives a well- 
defined recursion procedure for computing all coefficients Unm in a- finite number 
of steps. 

As mentioned, the method described above is a somewhat simpler variant of 
Sutherland's original method [Sut72] who, instead of the monomials M„, used 
somewhat more complicated basis functions which we now describe. Let 

L"/2J , 

(16) Hr.i.)=E(-^ri^^s^^r-^' 

denote the Hermite polynomial of order n € Nq, satisfying the differential equation 

i-dl + 2xd,)Hn{x) = 2nHn{x). 

( /Tl 

Let Mn denote the symmetric polynomial 

M^\yL) ^ Mn{Hr,^{xi),.. .,Hnr,{xN)) ^ ^ iJ„i (xp(i) ) . . . i7„„ (xp(A,) ) . 

These symmetric polynomials are obviously eigenstates of the differential operator 
H in for the free case A = with eigenvalues 2|n|. The key identity is [Sut72] 

^—(9, - dy) iHnix)Hm{y) + Hn{y)Hra{x)) 
X — y 

(17) n m 

= ^^C,.s(n,m) {Hn-r{x)Hm-s{y) + Hn-r{y)Hrn-s{x)) 
r—1 s—1 

for all n,TO € No and certain real coefficients Crs- Since a proof of this iden- 
tity is not contained in Sutherland's paper [Sut72] we provide a sketch thereof in 
Appendix ^ This identity shows that the action of H on the symmetric poly- 
nomials is triangular, which suggests that there are eigenfunctions — 

Mn^^ + J2in<n ^nmAfm^'' of H with eigenvalues 2|n| and a recursive procedure to 
compute all coefficients Wnm from Vnn = 1, as above. However, the explicit formulas 
for the coefficients Crs{n,m) were not provided in [Sut72], and they indeed seem 
rather difficult to obtain: we neither found them in the literature, nor where we 
able to derive them. 

It is interesting to note that the recursion relations in H15|) can be inverted to 
yield explicit formulas for the coefficients Wnm: introducing = 1 and a linear 
operator R acting on these coefficients as follows: 



^^nkOkm 

m<k<n 



we can rewrite the recursion relations as 

It follows that they can be inverted according to 



w„m = (l-i?)-^(5„,m-^i?^<5„ 



s=0 

where it is important to note that the expansion of the geometric series is well- 
defined since it only contains a finite number of non-zero terms; see below. From 
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the definition of tlic linear operator R given above and the fact that — Em = 
2(|n| — |m|) now follows that 

—A 1 - t^n.m I \^ &nki&kik2 ■ ■ ■ ^k^m \ 

U»j Unm - f)„,m + ^| _ |^ I ^nm + 2^ 2^ j _ jj^ I • 

^' \ s=l m<ki<---<k,<n llr=lVI I I r\J/ 

The restrictions imposed by the inequality in the second sum clearly implies that 
this series representation for the coefficients Unm only contains a finite number of 
non-zero terms. Also note that each term is well-defined. However, this formula 
is not very useful since the 6nm are not given by a simple formula. Indeed, to 
deduce these latter coefficients from Ijlbf) it is important to note that n — 2ej is not, 
in general, a partition, e.g. (3,2,2) — 2ei = (1,2,2). A similar remark appHes to 
the last term in (|13|l . We therefore implicitly used an extension of the monomials 
Mji to non-partitions. To make this precise we introduce an ordering symbol as 
follows: for each a G Nq we let p[a] denote the corresponding partition obtained 
by permuting the elements of a, e.g., p[(3, 1, 0, 4)] = (4, 3, 1, 0). We can then define 
Ma :— iWp[a]i which naturally extends the definition of the monomials Mn to non- 
partitions. Using this definition we deduce from (|13l) that 

JV L^^^J 

bnm = - '^n.jinj - l)(Sp[n-2e,].m ^ ^ X! 

j—1 j<k i^— 

X (2 - 52u,nj-nk) ((1 — ^i/,o)"-j " ^k) ^p[n- (!y+l)ej l)efc] ,m- 

Inserting this in (|18|) one hopes that, due to the Kronecker deltas, the sums sim- 
plify considerably. However, the appearance of the ordering symbol p[-] makes the 
resulting formula awkward to use. We therefore conclude that the Sutherland al- 
gorithm does not lead to simple explicit formulas for the eigenf unctions. The same 
difficulty arises in Sutherland's original algorithm described above. 

3. Alternative solution algorithm 

In this section we present our alternative method for solving the Calogero model 
defined by the Hamiltonian in and, in particular, prove Theorem ll.il 

We will to a large extent use the notation introduced in the beginning of Sec- 
tion|21 with the important difference that elements m, n e now will be ordered 
as follows: 

JV JV 

m ^ n ^ mk < ^ rik, yj — 1, . . . , N. 

k=N+l-j k=N+l-j 

3.1. A remarkable identity. We start by proving a particular functional identity, 
which is the starting point for our construction. 

2(1 - X)N and 



Lemma 3.1. Let cn 

Then 



(19) i?(x)F(x, y) = [H{y) + cn]F{x, y), 

where H = i?(x) is the Hamiltonian in 0) and similarly for H{y). 
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j=l ^ ■> ' j<k ^ ^ 



Proof. We set A/" = 2N , Xj = Xj, Xn+j — Uj, rrij — +1 and niN+j = —1 for 
j = 1, 2, . . . , TV. Then i7(x) - H{y) = n{X) with 

(20) n^Y. + + E ^ "^•^^^"^"'^ 

and F{x,y) = ^'o(X) with 

AT 

(21) *o(X) = Yl e-^""'^? Y[{Xk - Xj)"'^""'^. 

3 = 1 j<k 

To prove the lemma will we show by explicit computation that 

(22) (H-fo)*o(X) =0 
with the constant 

(23) fo = AK]mJ +Y,{l-\m]) 
and note that So = cn- For that we introduce the operator 

j=i J 

with 



[Xj - Xk) 



Note that ^'o = for all j = 1 , . . . , A/", and hence that 

V^a = 0. 

This implies the identity in H22|l since V = Ti — 8q, as can be shown by straightfor- 
ward computations. Indeed, 



N 



where the reminder terms 

N 



l + 2^Am,-mfc J + ^ ' ' ' 



j=i ^J-^k {Xk - Xj){X, ~ Xi) 

add up to the constant £q: upon symmetrization the double sum becomes indepen- 
dent of the Xj and equal to J2j^k ^t^jnik = A(^^ ^ ^ X]j ^^"^ the triple 
sum vanishes, as can be seen by symmetrizing in the indices j,k,£ and using the 
identity 



{Xk — Xj){Xj — Xe) {Xk — Xj){Xe — Xk) {Xj — Xe){Xe — Xk) 

□ 
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Remark 3.1. It is easy to see that we have, in fact, proved a more general result: 
the identity in (|22|l holds true for all A/" = 2, 3, . . . and arbitrary real parameters 
rrij. Obviously, one particular consequence of this latter result is the fact that the 
function ipo in ^ is the groundstate of the Hamiltonian iJ in Q . It is not difficult 
to see that we have proved a similar fact for a more general case: if all rrij > then 
Ti in ^20\) defines a self-adjoint operator on the Hilbert space L^{R^) with^oi^) 
l^21\) as groundstate. This is true since then is the Hilbert space adjoint of QJ, 
implying that T) = Ti — £q \s a. sum of non-negative terms {Qj)*Qj /rrij and thus 
defines a unique non-negative self-adjoint operator via the Friedrichs extension; see 
e.g. [RS75]. We thus recover a known generalization of the Calogero model where 
the particles can have different masses rrij > and such that its exact groundstate 
and groundstate energy can be computed exactly [For92,MMS03] . Other interesting 
special cases will be discussed in Remark 15. 21 

3.2. Integral transformation. The idea is now to apply to the identity in 119|l 
an integral transform II^i (/c '^vA'^^^Vj)^^ 4'nj{Vj)^ with the integration paths 
in |SJl and certain functions 4>j{yj) to be chosen such that this transform is well- 
defined. We observe that 



which shows that if we choose 



.,2 r) ^ 



(24) cjj.iy,) = e-^y.y;\ fi, = n, +\{N + 1~ j) 

with integers nj, then this transformation is well-defined for all A > 0. Indeed, for 
all n e Z^, 

^n(x) := n ^^"""^'^/j"') ^(^'y) = V'o(x)/„(x) 

with /„ the functions defined in @. 

The application of this integral transform to the l.h.s. of the identity in H19() 
obviously gives HF^lx.). To compute the integral transform of the r.h.s. we observe 
that 



(-C + y|)e-^^?yr"' - e-'^y' ((2n, - l)yf - (n, - l)(n, ~ 2)^;^-^) 

OO 



and that 



{yj-ykf yl{l-yj/yk)' 
for all \yj \ < Using these two facts and the shorthand notation 

7 = 2A(A- 1) 

we obtain by straightforward computations that 

N oo 

(25) HFn - EnF„ - ^(nj - l)(nj - 2)F„_2e, +7^111 '^4-(l-.)e,-(l+.)e,, 
j = l j<^kv=l 

where we used that X^jLi 2(2rij — 1) + cn — En, as given in 
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3.3. Construction of eigenfunctions. Equation (|25|l sliows that tlic action of 
H on tlie functions fn has triangular structure: HF^ is a linear combination of 
functions -Fm with m ^ n. Similarly as in the Sutherland algorithm, this suggests 
that the Calogero model has eigenfunctions of the form 

(26) Vn = a„(n)F„ + ^ a„(m)F^ 

with eigenvalues i?n and certain coefficients Q!n(m). Indeed, inserting this formula 
for V'n in H25I) . we obtain by straightforward computations that 

/ N 
Hilary = ^ E^an{m) - ^(mj + l)mja„(m + 2e^)+ 

m^n \ J — 1 

oo \ 

7^^^an(m+ (1 - :/)ej + (1 + j/)efe) jF^n- 

We conclude that the validity of the Schrodinger equation H^-^ — E^ip^ is implied 
by the recursion relations 

oo 

(27) 2(|n| - |m|) a„(m) = ^ ^ 9jk{'^; m)«„(m + E^^,), 

j<k 1^=0 

with gjk{i';m) and Ej^ defined in Q and ©, respectively; we used En — E^i — 
2(|n| — |m|). We now construct an explicit solution of (|27l) . The triangular structure 
of the eigenfunctions implies that Q;„(m) = unless m -< n or m = n, and that we 
can set 

an(n) = 1 

without loss of generality. This implies that the recursion relations in 127|) can be 
written as follows: 

an = <^n + San, 
where the operator S is defined by 

(28) (5a„)(m) ^ ^ 5^.,(^; m)a„(m + E^^,) 

^' ' ' ^' 'j<k v=0 

for m ^ n, which allows us to suppress the common argument m. This later 
equation can now be solved to yield 

s=0 

where the latter expansion of the geometric series is well-defined since it only con- 
tains a finite number of non-zero terms, as shown below. Using we deduce 
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that 

fqs^ vrr,^ _ \^ V" 9jskA''s;ni) ^ ^ g3.-ifc.-i(i^s-i;m + Ej;feJ 

(A„dJ(m) - 2. 2. 2(|n|-|m|) ^ ^ 2(|n|-|m + E^^ |) 

TT ^irfcrC'^r; n — ^^^^ E ) 

= Z.---Z. 2. ii 2(|n|-|n-V'- E-- n '^^(m + E.^.E t^^) 
ji<fci j,<fc, 1.1,...,^, r=i "^"^l"! I" 2^e=i^jtkt\) 

for all m ^ n. We now observe that 

2(|n|-|n-ELiE;;,J) = 4r 

and thus obtain {Tj). 

3.4. Properties of the reduced eigenfunctions. There remains to prove that 
the reduced eigenfunctions in Theorem 11.11 indeed are well-defined symmetric 
polynomials. We do this in three steps: we first establish that the functions /„ 
are symmetric polynomials, then, that the P„ are finite linear combinations of the 
functions /„, and finally, that all the expansion coefficients an(ni) are well-defined. 

A proof of the first fact can be found in [LanOl], but for the convenience of the 
reader we give the complete argument. 

Lemma 3.2. The functions /n are homogeneous symmetric polynomials of degree 
|n| and non-zero only if 

nj + ... + nN>0, Vj = 1,. ..,iV- 1. 

For each n G and each partition m, let 



where the second sum extends over all non-negative integers Kij and Vrs such that 

N j-1 N N 

1=1 1=1 1=1 

Then 

(29) fn= Yl PnmM^. 

m| = |n| 

Proof. Since \yj\ < \yk\ and \xj \ < \yk\ along the integration paths in /„, the terms 
in the fraction contained in its integral kernel can be expanded in binomial series 
in yj/yk and Xj/yk, respectively. The integrals can then be computed using the 
residue theorem, and this yields 

/n(x)=E(-l)— C')^' 
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where the last sum is to be taken over aU non-negative integers Kij and Vrs such 
that 

j-l N N 

(30) rij - ^ Kij + ^ Kji - ^ uij = 0. 

;=i i=j+i 1=1 

Recalhng the definition of the monomials Mm we deduce (|29ll . To prove that the 
functions /„ are homogeneous of degree |n| note that the degree of each monomial 
A/m in the decomposition in (^5)1 is given by 

N N 

i,i=i j=i 

It remains only to prove that the functions /n are polynomials, i.e., that H^Ofl only 
has a finite number of solutions. To this end consider the equation for j — N , 

JV-l N 

riAT = ^ Kin + ^ vm- 
1=1 1=1 

It is clear that it only has a finite number of solutions for each fixed njv- Now 
observe that, for j < N, the possible values of kji, I > j + 1, are determined by the 
equations with larger values of j. Also observe that the equation for each fixed set 
of Kji, I > j + 1, and rij has only a finite number solutions. The statement thus 
follows by induction in j, starting with j = N. □ 

Observe that either the symmetric polynomial /m or the coefficient an(m) is zero 
unless rrij + ■■■ itln > ior all j — 1, N — 1 and m ^ n. Clearly only a finite 
number of m € fulfill these conditions. This proves the second fact: the reduced 
eigenfunctions are finite linear combinations of the symmetric polynomials /„. Also 
note that all sums in the explicit representation in ([TJ) of the coefficients a„(m) 
truncate after a finite number of terms, and that they therefore arc finite. It follows 
that the reduced eigenfunctions indeed are well-defined symmetric polynomials. 
This concludes the proof of Theorem ll.il 

3.5. The one-particle case. It is interesting to note that Theorem 11.11 is non- 
trivial already in the simplest case N = 1. Since the Calogero model for iV = 1 
reduces to the harmonic oscillator with well-known eigenfunctions given by the 
Hermite polynomials Hn (see e.g. 22.6.20 in [AS65]), Theorem 11.11 in this case 
implies that the functions P„ in (|14(l are equal to the Hermite polynomials up to 
normalization. Comparing with the standard definition of the Hermite polynomials 
in (|16|) we obtain that 



l^Jn „_n ^ *■ J y > X 



' -|y|>|:r| 27ri?/ (1-x/y)^' 

where (z)„ denotes the Pochhammer symbol 

(z)o = 1, {z)n^ z{z + l)...{z + n~l), 

defined for z € C and n G Nq. Note that the series above truncates after a finite 
number of terms and thus is well-defined. This identity has an interesting inter- 
pretation. Observe that the Hermite polynomials can be generalized to arbitrary 
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complex parameters a as follows: 

It is straightforward to verify that this series reduces to a Hermite polynomial when 
a is a non-negative integer, and that it satisfies the Hermite differential equation 
{d^ — 2xdx + 2a)H^{x) = in the sense of formal Laurent series. However, it is 
important to note that the series defining Ha (x) does not converge anywhere in the 
complex plane but is only asymptotic unless a is a non-negative integer. Using this 
formal Laurent series we can formally rewrite the above identity as follows: 

TT / \ / dy 1 

2^ J\y\>\x\ 27ri (y-Xy 

for any complex A. For non-integer A the r.h.s. can be made well-defined by ex- 
changing the order of integration and summation. For integer X = m + 1 > 2 we 
can use the residue theorem to compute the integrals and recover the well-known 
identity 

Hn{x) = —, ^-r—Hn+mix) 

2™ (to + n)! ax™ 

obeyed by the Hermite polynomials, and we therefore obtained an interesting gen- 
eralization of this to the cases when n is not a non-negative integer. The inte- 
gral transforms in our identity looks like a fractional integral transform which, as 
is well-known, shift parameters of hypergeometric functions; see e.g. Chapter 13 
in [EWT53]. However, the details of our identity seem different. 

3.6. Alternative formulas for the eigenfunctions. The results above can now 
be used to construct another explicit series representation for the eigenfunctions 
of the A^-body Calogero model. For that we find it convenient to use a somewhat 
different normalization for the formal Laurent series Ha, 

(31) Pa{x) = 2-'^Ha{x)=f2c,{a)x''-'^, c,{a) = {-iy^-^^^±^. 

s=0 

The idea is to apply a particular integral transform to the identity in Ijiyi) which 
differs from the one in Section lX^ in that the simple powers yj' in 124|) are replaced 
by the formal Laurent series {yj ) . This leads to recursion relations which are 
somewhat different from those in the proceeding discussion but also can be solved 
explicitly. The advantage is that this recursion becomes trivial in the free case 
A = but, as we will see, it becomes somewhat more complicated to deduce. 

To obtain the recursion relations we need an explicit formula for x"pa{x), n e Z, 
as a linear combination of pa'{x), a' < a + n. 

Lemma 3.3. Let n G Z and a E C Then 

oo 

(32) x"-pa {x) = 5s (n, a)pa+n-2s {x) 



s=0 



with 



bo{n,a) = co(a) = 1, 



s-l 

bs{n, a) — Cs(a) — bj{n, a)cs^j{a + n — 2j), s > 0. 

3=0 
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and Cs{a) as defined in 

Proof. Observe that the definition of pa imphes that 



s=0 

in the sense of formal Laurent series. It follows that 



-2s 



x"pa{x) = co{a)pa+n{x) + ^ (cs(a) - CQ{a)cs{a + n)) x^^"" 

s=l 

oo 

= bo{n, a)pa+n{x) + ^ (cs(a) - bo{n, a)cs{a + n)) a;°+"" 



2s 



The statement now follows by repeating this procedure of breaking off the leading 
term. □ 

Remark 3.2. Although it is not evident from the statement and proof of Lemma lH.HI 
the series in ^XZl truncates for non- negative integers n at s = n, i.e., 

bs(n, a) = 0, s > n, n G Nq. 

This can be deduced by observing that differential equation solved by Pa implies 
the three term recursion relation 

2pa+l{x) - 2xpa(x) + apa^l{x) = 0. 

From this the truncation of the series in H32|l for non-negative integers n follows. 
This also shows that the series does not truncate for negative integers n. 

Explicit formulas for the coefficients bs{n, a) can now be obtained by solving the 
recursion relations in Lemma 13.31 

Corollary 3.1. Let s e No, n e Z and a G C. Then 
,(a- 2s+ l)2s 



bs{n,a) = i-iy 



I Vf_iV+^ V (« - + l)2s, (g + n - 2s + l)2(s-s,) 
+ 2^^ > 2^ 4''s,,!(s,_i -s,)!---(s-si)! 

j=l 0<Sj<---<s J V J i J' \ ^1 

Proof. Lemma l3 . 31 implies that 

s s-1 si-1 Sj_i-1 

6s(n,a) = Cs(a)+^(-ir ^ ^ c,^ (a) 

j — l Si— S2 — Sj— 

Csj-i-sj {a + n- 2sj) . . . Cs^si [a + n ~ 2si). 

The statement is now obtained by using the explicit form of the coefficients Cs(a) 
and simple properties of the Pochhammer symbol. □ 

Suppose that yjv > Vn-i > ■■■ > TJi > rnaxfc(a;fc). It is then clear from the 
previous discussion that the product 



j=i n^fc=i(i-2;/c/yj) 



X 
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is a formal multi- variable Laurent series in the variables yj. Appealing to the 
residue theorem we define 



(34) /r(x):=nf/ 



dyj , . A(7-A'-l) \ nj<fe(l ~ 2/j72/A:) 



A 



as the coefficient of the term {yi . . .yp^)^^ in this Laurent series; we use the super- 
script '{H)' to distinguish these from the analog functions defined in Section IX^ 
This prescription amounts to interchanging the integrations and summations. Fol- 
lowing the proof of Lemma 13.21 it is readily verified that the are well-defined 
symmetric polynomials. Similarly as in Section l3.3l we now obtain that 



oo oc 



j<Ck I'—l r.s—0 

X F^"^ 

n-(l+2r-i/)ej-(l+2s+iy)efc 

for Fn (x) = ■(/;o(x)/n (x); we used the fact that hp^, — 2vpi,, the functional 

identity in H19|l . and Lemma [3.31 As before we conclude that the action of H on 
" ( H^ 

the functions F^ has triangular structure, and that there are eigenfunctions of 
the Calogero model which are of the form = /9n(n)FA + X^m^n ■ 
The Schrodinger equation Htpn — Enipn is now implied by the recursion relations 

with 

^ oo 

(5W/3„)(m) = — — — 5^ ^ g*]J(:.;m)/3„(m + E*r), 

3<k u.t.u—ij 

where we have introduced the notation 

(35) g^^{v] m) = 2A(A - Vjvhtiv - 1; -t- 1 2t - v)hu(,-l - v;mk + I + 2u + v) 
and 

(36) ^^'^ = (1 + 2t - v)bj + {1 + 2u + v)ek. 
Using 

2(|n| - In - ELi I) = 4 ^(1 + + u,), 

1=1 

computations similar to the ones in Section lead to the following result: 
Theorem 3.1. For n e let 

P„=^/3„(m)/^f) 

m 

with the functions fm^ defined in and 

oo oo oo 

/3„(m)=5„(m)+$: E •■•E E 

S = l ji<ki tl, Ml, 2^1=0 is<fcs ts, Ms, 1^3=0 



4ELi(i + ^i + "0 
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with the quantities given in 1.^,5)) . ^.¥6)) and iV'jj) above. Then Pn is a reduced poly- 
nomial eigenfunction of the Calogero model corresponding to the eigenvalue Ea in 



Remark 3.3. The polynomials Pn constructed here should be identical with the 
ones obtained in Theorem ll.il We have checked this for N = 2 and various cases 
n e using MATHEMATICA. 

4. Solutions for the remaining classical root systems 

In this section we show that the constructions of the previous section can be 
adapted to the Bjf variant of the Calogero model defined by the Hamiltonian in 
()10|l . The construction is very similar to the A^v-i case, and we therefore are rather 
sketchy and concentrate on the necessary changes. To simplify notation we denote 
corresponding quantities in the Ajv_i and Sat cases by the same symbol, e.g. ■00, 
P(x, y), /„ etc. have a different meaning here and in Section|31 However, since the 
parameter fi will play a special role, we will write the Pjv Hamiltonian in (|10|l as 
H^. Moreover, the analog of the Hermite polynomials are the Laguerre polynomials 
denoted by the usual symbol Lii\x); see e.g. [AS65]. 

4.1. A remarkable identity. The analog of the key identity in Lemma [3. II is 



It is important to note that this identity now involves two Hamiltonians with dif- 
ferent coupling parameters fj, and A — /i. 

The proof is similar to the one of Lemma |3.1I by direct computation one can 
check that 




(39) 




j=i 



and the differential operator 




fijiHj - 1) 




+ H (X^-X^)^ (2(Am,mfc - 1) 




obey the identity in (|22|l with 



(40) 




The identity to cancel the three-body terms is now 




0. 



Equation (|S7|) is obtained from this as a special before. 
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This also gives, as another important special case, the groundstate 

N 

(41) = n (e"^'^'^') lii^i - ^')' 

j=l j<k 

and corresponding groundstate energy Eq of the -Bat Calogero model, which are of 
course well-known; see e.g. [OP83]. We note in passing that we have obtained a 
generalization of the Calogero model to particles with different masses, together 
with its exact groundstate and groundstate energy. 

4.2. Integral transformation. The kernel of the integral transform is now taken 
to be 

[[ 

with the same fij as in H24|l . Note that we need to restrict to even integers 2nj to 
get non-zero results; see (1^ below. 

Applying the resulting integral transform to the identity in (|37|l straightforward 
computations lead to 

N 

Hf^Fn = EnPn - ^ 2(2(n, + M - A) - - l)/'„-e, 

oo 

+ 4A(A - 1) ^ ^(2^. - l)F„-(l-.)e,-.e, 
j<k v=l 

with 

(42) £;„ = 4(ni -1-712 + ---+njv)+i^o, £^0 = A^(l + 2^ + 2A(iV - 1)) 



and = ipofn, 




_dyj ^2n,\ Ilj<ki^-yj/yh 



Note that the function /„ by a change of variables, x'j — *■ Xj and y| — *■ yj for all 
j ~ 1, . . . , N , becomes identical with the function /„ used in the An-i case. Thus 
Lemma H-j. 21 directlv extends to the Bn case. 

4.3. Construction of eigenfunctions. The action of the Bn Hamiltonian on 
the functions Fn is triangular, which suggests that the Hamiltonian in (|10|l has 
eigenfunctions of the same form as in H26() with eigenvalues in (|42|l . As before, 
the Schrodinger equation Hijjn = E-a'ipii is implied by the recursion relations in (|27|l . 
where now 

93k{v- m) - 4A(A - 1)(2;. - 1)(1 - Vo)(l - 5,k) 

(44) 

- 2{2{rhj + fi- X) + l)rhj6^flSjk 

and 

(45) E^-, -(l-z/)e,+z.efc. 

With these substitutions the arguments in SectionOlgo through as they stand, and 
we obtain the analog of Theorem II. II for the Bn case. To summarize, the reduced 
polynomial eigenfunctions of the Bn Calogero Hamiltonian in are given by the 
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functions Pn in Theorem M.lX where J^, |^ and have to be replaced by {4-S\j , 
\4-4\j and i45}) , respectively. The corresponding eigenvalues are given by i4S\ ). 

4.4. The one-particle case. The Calogero Haniiltonian in (|10|l for = 1 

reduces to —d^ + x'^ + — l)x~^ with well-known eigenfunctions given by the 
Laguerre polynomials with non-degenerate eigenvalues; see e.g. 22.8.18 in [AS65]. 
Comparing this with our eigenfunctions we can conclude that, in the case iV = 1, the 
polynomials P„ (x) constructed above are proportional to the Laguerre polynomials 
^(^^-i/2)^_^2^^ By straightforward computations we obtain that 

r (M-i/2) _ „i-(i-A) T{n + X) r dy^ , . I 

(A)„ J\y\>\x\ ^TTi ^ [y - x)^ 

where we have made the substitutions y^ ^ y and — > x as well as extended the 
definition of the Laguerre polynomials to arbitrary complex values a as follows: 

T (m-1/2) _ V cJ-^i'^-'^) (g -f M + 1/2 - fc)fc a-k 

r{a + l-k)k\ 

with the Gamma function F; note that La^ ^^'^\x) is a formal Laurent series obey- 
ing the Laguerre differential equation {xd^ + {^+l/2—x)dx+a)Li'^~^^'^\x) = 0, and 
it reduces to a Laguerre polynomial when a is a non-negative integer; see e.g. 22.3.9 
in [AS65]. As in Section f3.5l the r.h.s. in the previous equation has to be interpreted 
by exchanging integrations and summations as well as computing the integrals of 
the individual terms using the residue theorem. For integer values A = m -f- 1 > 2 
we recover the well-known classical identities 

Lir^/^Hx)^i~ir^Lirl^'-'-H-)- 

As in Section [3. 51 the general case is similar to known identities involving fractional 
integral transforms; see e.g. Chapter 13 and 16.6.(5) in [EWT53]. 

Using these results it is straightforward to extend our alternative formulas for 
the eigenfunctions in Section rTHl to the Bn Calogero model. 

5. Concluding remarks 

In the present paper we extended a solution method for the Sutherland model 
[LanOl, Lan05] to the A^v-i and B^ Calogero models. Below we discuss various 
open questions and mention interesting complimentary results. In particular, in 
Remark 15.11 we comment on the issue of whether our solution is complete or not, 
and in Remark 15.21 we point out further identities generalizing our key results in 
(I19|l and H37|l . In the concluding Remark 15.31 we sketch an interesting alternative 
interpretation of our method. 

Remark 5.1. The main difference between Sutherland's method [Sut72] and ours 
is that he expands the reduced eigenfunctions of the Calogero model in the mono- 
mials Mn defined in H12|l . whereas we obtain the reduced eigenfunctions as linear 
combinations of the more complicated functions /n in 0. 

The monomials Mn labeled by partitions n = (rj-i, n2, ■ ■ ■ ^un) are a basis of the 
space of symmetric polynomials, and it is therefore obvious that the eigenfunctions 
obtained by Sutherland's method are complete. 

On the other hand, the functions /„ and also our eigenfunctions are labeled by 
unrestricted A'^-tuples n e Z^, and it therefore seems that we are working with 
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overcomplete sets of functions. While this is the reason why we can get more 
expUcit formulas, it also makes the questions of completeness of our solution more 
complicated. It would therefore be interesting to extend our results to the higher 
differential operators commuting with the Hamiltonians since this might shed some 
light on this important issue. 

It is interesting to note that the very same functions /„ appeared as building 
blocks in all the models we so far have solved by our method: in the Aj^^i Calogero 
model they appear as functions of the variables Xj (see (O), in the Bn Calogero 
model as functions of a;| (see (031), and for the trigonometric Sutherland case we 
found the very same functions but in the variables zj = e'^^ [LanOl]. 

We finally mention that we convinced ourselves that, for A = 1, the are up to 
a difference in sign identical with the Schur polynomials. 

Remark 5.2. The key to our method was the identity in Lemma l3.1l but in its proof 
we obtained a more general result, given in H2()|l - H23|l . which has other interesting 
special cases, as discussed in Remark [3.11 We now point out further interesting 
identities which can be obtained from this general result. 

For example, it is possible to generalize Lemma 13.11 by allowing the particle 
numbers N and M in the x- and y-variables to be different: choosing Af = N + M, 
Xj ~ Xj, rrij = +1 and Xpf+k = Uk, ruN+k = —1 for j = 1,2,...,N and k — 
1, 2, . . . , M we obtain the identity 

(46) ffjv(x)FAr,M(x, y) [i?A./ (y) + C7v,m]^jv,m(x, y) 

for the function 
(47) 

1^2 



f \ -ix^ TT lv'^,'Wl<j<k<Mi^k Xj)'^Y\^^j^j^^j^j{yK yjY 

^A,M(x,y) = J_J_e 2^,||e2y./ 



j=i Y{j=iY{K=i{yK - x-jY 

oi N + M variables and the constant 

(48) CN,M = KN-Mf + {N + M){l-X)- 

Hn{'k) in H46|) is the Calogero Hamiltonian in ^ and similarly for iJAf (y), where we 
now also have to indicate the particle numbers and variables. Another interesting 
family of identities is obtained by choosing Af = N + M , Xj = Xj, rrij = +1 and 
Xiq+k = Uk, niN+k = 1/A for j = 1, 2, . . . , and k = 1,2,..., M. This yields the 
identity 

i?jv(x)FAr,j\/(x,y) = [-AiJAf4/A(y) + CN,M]FN,M{^,y) 

for the function 

N M 

1^2 -n- _ 1, ,2 , 



(49) 



i^A,Af(x,y) = n e^^'^' n e"^"''^^ n - ^j-)' 

j=l J=l l<j<k<M 

N N 



1<J<K<M 3 = 1 K=l 



and the constant 



(50) CNM = (A^A + M/\f + A^(l - A) + M{1 - 1/A); 
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the Calogero Hamiltonian for the variables y is now given by 

^ 1 

(51) i?M4/A(y) = E + yj/>^') + 2(1/A)((1/A) -1)Y^ ^ 

j=i j<K ^y-' 

We thus recover a well-known duality between Calogero models with reciprocal 
coupling constants. 

It is interesting to note that corresponding identities also exist in the Bpf case, 
but we do not write them down since they can be obtained from our general result 
quoted in H,'W|I -14() () in the same way as explained in the An-i case: to obtain the 
analog of 147|l - l|48() fix the parameters as in the An^i case and, in addition, set 
fij = /i, fiN+k = A — /i, and for the analog of (|49(l - (|51|l the additional parameters 
are to be fixed as fXj = /z, fJ-N+k — — 1). 

With these identities one can obtain many more explicit formulas for the eigen- 
functions of the Calogero models using our method. It is interesting to note that 
similar identities were previously found also in the Sutherland model and its elliptic 
generalization using quantum field theory techniques [Lan04], and we believe that 
such identities should exist also for other integrable many-body systems. 

Remark 5.3. We now sketch an interesting alternative interpretation of our method 
which, when explored in more detail, could shed more light on the questions dis- 
cussed in Remark 15.11 ^ Suppose that we want to construct eigenfunctions of the 
Calogero Hamiltonian in ^ of the following form 

N 

Xn = XoPn, Xo(x) = Yl e"^"'^' 
i=i 

with Pji linear combinations of arbitrary monomials 

„n — ^"1 ^"2 . . . 

-A. . N ' 

where the fij — rij + Sj are integers shifted by common amounts Sj which, at this 
point, are arbitrary. If we ignore the issue of square integrability it is easy to 
construct such eigenfunctions: with the ansatz 



P„(x) =x"+ J2 «n(m)x" 



we find by straightforward computations that the Schrodinger equation Hxn = 
EnXn is equivalent to = |n| + |s|, |s| = si -|- S2 -I- . . . -I- sat, and the recursion 
relations in H27() for the coefficients an(m). As shown in Section lX^ these relations 
can be easily solved, and their explicit solution is given in Theorem ll.il It is impor- 
tant to note that the construction of these eigenfunctions Xn is not restricted to the 
Calogero model, but one can easily generalize it to construct similar eigenfunctions 
for non-integrable systems like the generalized Calogero model where the particles 
have different masses. However, these eigenfunctions are completely uninteresting 
from a physical point of view: the series in the previous formula do not converge and 
are only asymptotic. The fact that makes the Calogero model special is that there 
exists an operator T which maps these unphysical eigenfunctions to physical ones, 



We would like to thank Vadim Kuznetsov for useful discussions on this point. 
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and this operator is given by the function i^(x, y) in Lemma |H.ll in the following 
sense: for an asymptotic series 

X(x) = Xo(x) a(m)x™ 

let 

^(x)(x) = 5]a(m) n ^) J^(x,y)xo(y)y* 

with the integration paths defined in ©. Note that this map is well-defined if we 
set Sj — X{N + 1 ~ j), and then the functions 

-^(^'nXo) = V'o/n 

are equal to the building blocks of our solution given in ^ and ((BJ. Moreover, 
since Lemma l3 . 1 1 implies that 

we conclude that 

V'n := ^iXn) 

is an eigenfunction of the Calogero Hamiltonian with eigenvalue — + cjv- 

Obviously the generalized identities pointed out in Remark IS . 21 provide operators 
which, similarly, transform unphysical eigenfunctions of the Af -variable Calogero 
model to physical iV- variable ones, and thus our results in this paper give the further 
explicit formulas for the eigenfunctions mentioned in Remark 15.21 

We note that the operator ^ is similar to the Q-operator which has appeared in 
a separation-of- variables approach to the Sutherland model [KMS03] . 
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Appendix A. The relation to Calogero's original model 

We give here a brief account of the relation between the Hamiltonian in and 
the following one studied by Calogero [Cal69,Cal71]: 

AT - 

Heal = - E < + E ( ^'(^^- - ^^)' + 2A(^ - 1)?^-^ 

It seems well known that they differ only in their center of mass motion^ but we 
have been unable to find a discussion of this in the literature. 
The simplest way to see this is to verify the following identity: 

Hcai + ^'^{xi +X2-\ h XNf = H for Nlj"^ = 1. 



'We thank F. Calogero for explaining this to us. 
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Since Hcai is translational invariant it is possible to use center of mass coordinates 
with ro = (xi + X2 + ■ ■ ■ + xn)/^/N and other coordinates rj , j — 1,2, . . . , N — 1, 
hnearly independent of tq and 

were Hcm only depends on the rj with j > 0. Wc thus conclude that 

where we used that Nuj'^ — 1. This makes manifest that the difference between 
the two Hamiltonians H and Hcai lies only in the center of mass motion: in the 
former it is trapped in a harmonic oscillator potential, and in the later it is free. 
Thus the Hamiltonian Hcai has discrete spectrum only if its center of mass motion 
is fixed, i.e., if only Hcm is considered. One the other hand, H has purely discrete 
spectrum even if the center of mass is not fixed, and it is therefore simpler to work 
with. 

It is instructive compute Hcm explicitly. For iV = 3 the center of mass coordi- 
nates are tq above and 

n = ~ ^2), r2 = "^(2^1 + 2^2 - 22:3). 

By a straightforward computation follows that 



Hcm = -dt, - 9^ + iuj\ri + r^) + 2A(A - 1) + 4 



2 



2rl [rl-irlY, 

This computation extends straightforwardly to the general case by a generalization 
of the center of mass coordinates to arbitrary N . The latter can be found in 
Reference [Cal71], for example. 

Appendix B. Proof of the identity in (fT7j) 

In this appendix we sketch a proof of the identity in l|17(l . It should be noted 
though that the proof does not give the explicit form of the coefficients Cr,s^ which 
seemingly is rather complicated. 

The key ingredient in the proof is the following: 

Lemma B.l. Let n, m e No and n > m. Then 

n — m— 1 



m— 1 _ n— l_m — 1 

-mx y — my x 



Proof. The case n = m is easily verified. Now suppose that n > m + 1. Then 



1 rf-m m 

-{d^ - dy){x'''y"' + = n— - 



X — y X — y 



x-y 

Expand the fractions in geometric series to obtain 



^m— 1 m— 1 
-m- ^ (a;"-™+i - 



n— m— 2 n — m 

En-2-k m+k _ n-l-fc m-l+k 

X y — m 2, 2; y . 

k=0 k=0 
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Collect terms of equal degree to deduce the statement. □ 

The next step is to use the series representation of the Hermite polynomials 
and apply Lemma [B. II to each term separately. Recall that, for each n € No, the 
Hermite polynomial iJ„ is a polynomial of degree n. Hence, this procedure clearly 
gives a polynomial of degree at most n — 1 in the variables x and y. The validity of 
the identity in IjlTfl now follows from the fact that as n runs through the integers 
less than or equal to m, for some m e Nq, the Hn form a basis for the space of 
polynomials with degree at most m. 
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